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(-H I Abstract. We study the parity violation in the cosmic microwave background (CMB) bispec- 

trum induced by primordial magnetic fields (PMFs). Deriving a general formula for the CMB 
bispectrum generated from not only non-helical but also helical PMFs, we find that helical 
PMFs produce characteristic signals, which disappear in parity-conserving cases, such as the 
intensity-intensity-intensity bispectra arising from X]^^^ = odd. For fast numerical calcu- 
lation of the CMB bispectrum, we reduce the one-loop formula to the tree-level one by using 
the so-called pole approximation. Then, we show that the magnetic anisotropic stress, which 
depends quadratically on non-helical and helical PMFs and acts as a source of the CMB fluc- 
tuation, produces the local-type non-Gaussianity. Comparing the CMB bispectra composed 
of the scalar and tensor modes with the noise spectra determined by the cosmic variance, 
we find that assuming the generation of the nearly scale-invariant non-helical and helical 
PMFs from the grand unification energy scale (lO^^GeV) to the electroweak one (lO^GeV), 
the intensity-intensity-intensity bispectrum for Yln=i = odd can be observed under the 

2/3 1/3 

condition that B^^^^B^!^^^ > 2.7 — 4.5nG with .BiMpc and i3iMpc being the non-helical and 
helical PMF strengths smoothed on 1 Mpc, respectively. 
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1 Introduction 

The cosmological parity violation is a key feature of ultra-violet completion of general rela- 
tivity and hence a lot of researchers have extracted their signals from several cosmological 
phenomena [1-5]. In particular, the effects on the cosmic microwave background (CMB) have 
been well-studied and the cosmological parity violation has been verified by analyzing the 
non-vanishing cross-correlated power spectra between the intensity (I) and i?-mode polariza- 
tion (B) anisotropies and those between £'-mode (E) and i?-mode polarization anisotropies 
[6-11]. Furthermore, beyond the linear-order eff'ects, the impacts of the parity violation on 
the graviton non-Gaussianities have recently been discussed [12-14]. According to Ref. [12], 
unlike the parity-conserving non-Gaussianity, the parity-violating one induces the signals 
arising from Yll=i 4 = odd in the GMB ///, HE, lEE, EEE, IBB and EBB bispectra 
and also those coming from Yln=i = even in the CMB IIB, lEB, EEB and BBB bispec- 
tra. In these correlations, the /// bispectrum from Yl\=i = odd is expected to bring in 
the detectable information of the parity violation. 

On the other hand, if there exists the primordial magnetic field (PMF), which is a fa- 
vored candidate for the seed field of microgauss-level magnetic fields in galaxies and cluster 
of galaxies [15-17], their power spectrum may involve the parity- violating component [18- 
22]. Like the above non-magnetic cases, this so-called helical PMF induces the characteristic 
signals in the CMB IB and EB correlations [23-27]. Although concrete limits on the mag- 
nitude of the helical PMF have not obtained yet, these studies imply that the IB and EB 
correlations are detectable if helical PMFs have nanogauss-level magnitudes (at the present 
time) and these spectra are nearly scale invariant. However, assuming the Gaussianity of the 
PMF, the beneficial signals are generated also in the CMB bispectra due to the quadratic 
dependence of the CMB fiuctuation on the PMF. In the case where only non-helical PMFs 
exist, the contributions of PMFs to the CMB bispectra have been deeply investigated in 
Refs. [28-36]. 

In this paper, we newly consider the effects of both non-helical and helical PMFs on the 
CMB bispectrum. Based on our computation approach [28-31], we derive a general formula 
for the CMB bispectrum induced by the non-Gaussianity of the PMF anisotropic stress 
coming from not only non-helical PMFs but also helical ones. Then, we confirm the existence 
of the foregoing parity-violating signals such as the /// bispectrum from Y^=i = odd. By 



- 1 - 



the pole approximation mentioned in Ref. [28], we reduce this formula to a form suitable for 
the fast calculation in the case where the non-helical and helical PMFs have the nearly scale- 
invariant spectra. In this process, it is shown that the bispectrum of the PMF anisotropic 
stresses has the local-type shape even if helical PMFs exist. Computing the CMB /// 
bispectra composed of the scalar and tensor modes and estimating the signal-to-noise ratio, 
we analyze how the helical PMF affects the CMB bispectrum and show how large the PMF 
strength is required for the detection of the III signals from X]n=i ~ odd. 

This paper is organized as follows. In the next section, we summarize the expressions and 
statistical properties of both non-helical and helical PMFs. In Section 3, the analytic formulae 
and numerical results of the CMB bispectra, and the signal-to-noise ratio are presented. The 
final section is devoted to the summary and discussion of this paper. Throughout this paper, 
we obey the definition of the Fourier transformation as 

and the rule that the subscripts and superscripts of the Greek characters and alphabets run 
from to 3 and from 1 to 3, respectively. 



2 Statistical properties of non-helical and helical magnetic fields 

Let us take into account the large-scale primordial magnetic field (PMF), i3^(x, r), which is 
generated in the very early Universe and behaves as a source of the CMB fluctuation, on the 
homogeneous background and small perturbative Universe as ds^ = a^[— dr^ + Ihobdrdx^ + 
{She + hi,c)dx''dx'^]. Here, a and r denote the scale factor and conformal time, respectively. 
Neglecting the effects of the back reaction of the fluid on the evolution of magnetic fields and 
considering the flux conservation, the PMF evolves as B^{x,t) tx l/a(r)^. Each component 
of the energy momentum tensor composed of the PMF is given by 



Svra^ 



i?2(x) , 



n (x'^) = n(x^) = , 

i?2(x) 



Tt(xn 



47ra4 



5^,-B^x)i?e(x) 



(2.1) 



The spatial parts in Fourier space are expressed as 



T'^iKr) ^ p,{t) [5'',AB(k) +n^,(k) 

AB(k) 



87rp^,oy (27r)3 



i?^(k')i?6(k-k') 



n 



1 



Bc\ 



4vrp^,0 J (2vr) 



d^k' 



i?^(k')Be(k - k') 



(2.2) 



where = p^^oo"^ is the photon energy density with p^^Q being its present value. After this, 
for simplicity of calculation, we ignore the trivial time-dependence and hence the index can 
be lowered by 5hc- 
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Conventionally, assuming the Gaussianity of the PMF, the power spectrum is expressed 
as [23] 



{Ba{k)Bb(k')) 



(2vr) 



PB{k)Pab{k) + iVabckPBik) 6{k + k') 



(2.3) 



where k is a unit vector, rjabc is the asymmetric tensor normalized by 7^012 = and Pafe(k) = 
^ab ~ kakb is a projection tensor coming from the divergence free nature of PMFs. The first 
and second terms in the bracket represent non-helical and helical contributions, respectively. 
For a mathematical relation: 

^lim^ (B(k) • B(k')> > ^lim^ k x B(k) • B(k'; 

the power spectra of non-helical and helical PMFs obey such a magnitude relation as 

Psik) > \PB{k)\ . (2.4) 

In order to formulate the CMB bispectrum, it is convenient to use a normalized divergenceless 
polarization vector in two circular states, e'^^\ as shown in Appendix A. Then, the above 
expression changes to 



{Baik)Bb{k')) 



(2vr)3 
2 

(2vr)3 



6{k + k') [PB{k)-aPts{k)]e^:Hk)elr''\k) 



a=±l 



6{k + k') ^ [Psik') - aPsik')] 6i"'^)(k')6r(k') . (2.5) 



(T=±l 



This implies that the second terms of the brackets in the first and second equalities creates 
the difference of the magnetic power spectra between two circular states as o" = ±1. 

To parametrize the magnetic field strengths, we introduce the quantities smoothed on 

r as 



B^ = (B(x) • B(x)) I, 
d^k f d^k 



(27r)3 J (27r) 



[Ba{k)Ba{k')) e--'('='+fc'')/2e^(k+k')-x ^ 



= r\ (B(x) • [V X B(x)]) 



d^k f d^k 



(27r)3 J (27r)3 
Assuming the simple power-law spectra as 

PB{k)=ABk''^ 

the spectral amplitudes are written as 

(27r)"s+5^2 



r\ik'kbVabc (i?a(k)i?e(kO> |e-^('^^+^^'^)/2e^(k+k')-x 



A 



B 



Peik) = Aek""^ , 

(27r)"B+5^2 



\Ab\ 



'nn+4 



+3 



(2.6) 
(2.7) 

(2.8) 

(2.9) 



where F is the Gamma function and kr = 27r/r. Note that unlike Ab, Aq can take both 
positive and negative values. Equation (2.4) leads to a constraint on the PMF strengths: 
\Br/Br\ < 0.288 if the PMF spectra have nearly scale invariant shapes as ub = = —2.9. 
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The PMF anisotropic stress, Usabi induces the CMB anisotropy and therefore we require 
their bispectrum for the computation of the CMB bispectrum. Considering Equation (2.5), 
this can be straightforwardly calculated as 



(nB„6(ki)nBcd(k2)nBe/(k3)) = (-4^p. 



■3 

.n=l o-n=±l 



X5(ki - k'l + k'3)5(k2 - k'2 + k'i)5(k3 - k'3 + k'2) 

+{a 6 or c o d or e o /}] , 



(2.10) 



where is the Alfven-wave damping length scale [37, 38] as A;^^ ~ C'(0.1)Mpc and the curly 
bracket denotes the symmetric seven terms under the permutations of indices: a f-)- 6, c f-)- d, 
or e ^ f . For the sake of avoiding the IR divergence, we limit the range of the PMF spectral 
indices as nB,nB > —3. 

Like the discussion in Ref. [28], if the tilts of the PMF spectra are enough red as 
~ ng ~ —3, the shape of the bispectrum of PMF anisotropic stresses depends strongly 
on the behaviors of the integrands at around three poles, namely, k'l^k^^k'^ ~ 0. In this limit, 
the bispectrum (2.10) reduces to 



{IiBab{^l)^Bcd{^2)^Bef{^3)) ~ (-4vrp^,0 



uAb ^^n^+3 87rl 
nB + 3 * 3 8 



E {PB{ki) - aMki)} {PB{k2) - <T2PB{k2)} 

a2,cT3=±l 



X5a,.6|,'^^^(ki)6(-'^^)(ki)6(-)(k2)6^^'^^^(k2) 

+ E i^Biki) - cTiPtsiki)} {PBiks) - a^PBiks)} 

(T1,(T2=±1 



xer^Hki)erf 



(ki)5,,,e(-2)(k3)e5.'^^^(k3) 



(-<^i)/ 

cTiPB(.k2)} {PBiks) - asPBiks)} 



+ E i^^(^2) 

(71,(T3=±1 

xe("'^^)(k2)e^)(k2)6r"^^(k3)6r^(k3)5c,/ 
+{a o fe or c O (i or e o /}] , 



(2.11) 



where k^, = lOMpc"'^ and a are a cutoff scale of the integrand and a parameter fixing the 
uncertainty of the amplitude associated with the approximation, respectively, and we have 
evaluated an integral at around each pole by following 



/ d'k' E {PB{k')-aPsik')}e(:\k')ei^^\k') 

J LI 



a=±l 



. (2.12) 

ub + 3 



In Equation (2.12), due to the summation over a = ±1, the contribution of the second term of 
the bracket vanishes. This implies that the effects of the helical PMF are tiny at around each 
pole. Prom Equation (2.11), we can see that the bispectrum of the PMF anisotropic stresses 
dominates at the squeezed limit such as /ci ~ /c2 J?> ^3 and has the identical /c-dependence to 
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the local-type bispectrum of curvature perturbations [39]. Thus, we conclude that with and 
without helical PMFs, the shape of the non-Gaussianity associated with the PMF anisotropic 
stress is classified into the local-type configuration. 

In the next section, we compute the CMB bispectra generated from the non-Gaussianity 
of the PMF anisotropic stresses. 



3 CMB bispectrum from non-helical and helical magnetic fields 

In this section, we investigate the effects of both non-helical (B) and helical (B) PMFs on 
the CMB bispectra. At first, based on the formalism presented in Refs. [28, 29], we derive 
their exact and optimal formulae generated from Equations (2.10) and (2.11), respectively. 
Next, through numerical computations, we analyze the magnitudes and shapes of the CMB 
bispectra and examine whether the parity-violating signals coming from helical PMFs can 
be detected or not. 



3.1 Formulation 

The CMB intensity and two linear polarization fields (X = I, E, B) generated from the 
scalar-, vector- and tensor-mode perturbations {Z = S, V, T) are expanded by the spherical 
harmonics as 



= Z^(^x'lm^i'M ^ (3-1) 



X 



where n is a line-of-sight direction. According to Refs. [40, 41], each spherical harmonic 
coefficient, a^x\m^ given by 



(3.2) 



where A = 0, ±1, ±2 expresses the helicity of the scalar-, vector- and tensor-mode perturba- 
tions, X = 0, 1 discriminates the parity-even (/, E) and -odd (B) fields, and ^(^) and Tj.^^{k) 
are the primordial perturbation and transfer function of each mode, respectively ^. 

If there exist large-scale PMFs, their anisotropic stresses generate additional fluctuations 
in the CMB. The PMF anisotropic stresses survive and become a source of the gravitational 
potential before neutrinos decouple and they are compensated by the neutrino anisotropic 
stresses. Thus, gravitational waves and curvature perturbations logarithmically grow even on 
superhorizon scales prior to neutrino decoupling and produce the CMB tensor- and scalar- 
mode anisotropics at recombination epoch [42-44] ^. In contrast, due to the decaying nature 
of the vector potential, the CMB vector-mode anisotropies are induced by not this mechanism 
but the vorticities of photons enhanced by the Lorentz force from the PMF anisotropic 



'^Here, we set O" 
■^Recently, the S' 

the inflationary and the radiation-dominated eras consistently [45, 46] 



■^Recently, the solution of the curvature perturbation is being reanalyzed by treating the effects in both 
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stresses [42, 47, 48]. Consequently, we can summarize the scalar-, vector- and tensor- mode 
initial perturbations as 



-Rj In 



TB 2 " 



(k)nBafe(k) , 

\i t> / ^ 

(^^)(k)nBa6(k) , 



(3.3) 



6i?^ In 



where O^^^ ,0^^^^ and O^^^^ are the projection tensors decomposing into the scalar-, vector- 
and tensor-mode variables, respectively, and defined in Appendix A. ^^^^ and correspond 
to the curvature perturbation and gravitational wave on superhorizon scales, respectively [42- 
44], and depend on the production time of the PMF, tb, the epoch of neutrino decoupling, 
Ty ~ lMeV~^, and the ratio between the energy densities of photons and all relativistic 
particles, i?^ ^ 0.6, for r < t^. As the upper and lower values of r^^, we take the energy 
scales of the grand unification (lO^^GeV) and electroweak symmetry breaking (lO^GeV), 



^^'^ and T^j , we use the standard 



corresponding to Tu/tb = 10^'' and 10^, respectively. As 7^ 
cosmological transfer functions independent of PMFs [40, 49-51] because the evolution of the 
cosmological perturbations are little-affected by PMFs posterior to neutrino decoupling [43]. 
On the other hand, as , we should use the form including the effects of the PMF on the 
cosmological perturbations shown in Refs. [42, 47, 48]. 

Using Equation (3.2), the CMB bispectrum is formulated as 



n 



(^n) 



An 



n ^LmS^ri 



(3.4) 



\n=l 



To obtain the explicit formula for the CMB bispectrum involving the dependence on PMFs, 
we have to compute the angular bispectrum of the initial perturbations, <^ 11^=1 ^l^m„(^"))- 
Expanding all angular dependence in the delta functions and the contractions of the projec- 
tion tensors and wave number vectors in <^nn=i ^e^rl (^n)^ ^P^^ spherical harmonics 
base on Appendix A, and expressing the angular integrals of their spherical harmonics with 
the Wigner symbols [29], this is obtained as 



\n=l 



ll 12 ^3 
mi 7712 ITT-S 



LL'L" ^,^',^"=±1 



ll 12 t3 

L' L" L 



n 

.n=l 



{PB{k'i)+SPB{k'i)} 



with 
eS"S\ 



E 

L1L2L3 



X {PB{k'2) + 5'Pb(4)} {Psik'^) + S"PB{k',)} 
^•/'l"L^i^(^3' ^i)/ll'^2^(^'i5 ^2)/f/f//^3^(fc2 5 ^3' ^s) 



y'^dyjis {r2,y)jL2 ir2y)jLi (ny) 



(3.5) 
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-L2+L3 ( 1 \ Q Q Q jQS"~S" tOS-S tOX~X 

^LiL2L3^LalL" ^ L2IL ^ L^12 



f - 2 (8vr)3/2ij^ In (t./tb) (A = 0) 



X < 



|(8J)3/2a (A = ±1) 

-4(8^)3/2i?^ In (t./tb) (A = ±2) 



Here, ji is the Bessel function and the / symbol is defined by 



JS1S2S3 



(2/i + l)(2;2 + l)(2/3 + l) / h h h 
47r V ■51 ^2 S3 




(3.6) 



(3.7) 



/ £2 i3 \ 

The confinement of mi , 1712 and m-i to the Wigner-S? symbol ( guarantees 

\mi m2 ms J ^ 

the rotational invariance of the CMB bispectrum. If = 0, this is consistent with the 
corresponding equation for the non- helical case [28-31]. Note that unlike Pb, Pb associates 
the spin as S,S',S" = ±1. As described in Equation (3.8), this differentiates the multipole 
configurations of non-helical terms from those of helical ones. Via the summations over 
S, S' , S" , Ai, A2 and A3, the CMB bispectrum from non-helical and helical PMFs is explicitly 
written as^ 

/ 3 \ 



\n=l 



(Zn) 



) CziC'zjC'zs (-4vr/3^o^ ^ 
nil 1^3 i ^ J V ;> 



E 

LL'L" 



h h h 

V L" L 




rO rO 7-000 
L 1 La is L'l L'/ L^' 





L' 


i2 ] 






L" 


^3] 








II 




^2 




[ 1 


1 








1 





27r2 



'X„lr 



(kn) 



k?dk[ [PBik[)Qil^^^, 



/"OO /"CXD /'OO 

/ A^dAjL.ikiA) / B^dBjL'ik^B) / dC j v^k^C) 
lo Jo Jo 

- PB{k'2)Q^;:lL,^^L] jL',{k'2B)jL'i{k'2C) 



PB{k'2)Q%r,r, 



„or01-l rOl-l rOl-l 7-01-1 7-01-1 7-01-1 
^°^L'^lL^L2lL^L'^ IL' ^ L'^ IL' ^ L3 IL" ^ L'^IL" 

y93-Ars7-0|^i|-|-^il7-0|-^2|-|A2|^0|A3|-|A3|^(e) „(e) „(e) 

^Li £1 2 t2 2 ^L'l ia 2 ^Li,ei,xi^L[/2,X2^L'{/3,X3 ' 



(3.8) 



Caution about a fact that |A| is determined by Z, namely, |A| = 0, 1,2 for Z — S,V,T, respectively. 
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where Ns is the number of the scalar modes constituting the CMB bispectrum and 
Cz= < 



i -^{87Tf/^R^ln{T,/TB) iZ = S) 

|(87r)3/2 {Z = V) ■ (3.9) 

-4(87r)3/2i?^ In (tJtb) {Z = T) 



Here, we have introduced the filter functions as 

2l^,L2,L = ('^^3,^+1 + ^L'^,\L~1\){5l2,L+1 + 5l2,\L-1\) + ^L's,Lh2,L , 

^lIl2,L = (^L'^,L+1 + ^L'^,\L-1\)^L2,L + SL'^,L{h2,L+l + (^L2,|L-1|) , (3-10) 

which come from the above summations and selection rules of the Wigner symbols [41] and 
ensure Lg + L2 = even, = odd and Li + £1 + xi = even, respectively. Considering these filter 
functions and a relation derived from the selection rules as 

3 

J2 i^n + Xn) +L2 + L'2 + 4' + ^3 + 4 + 4' = even , (3.11) 

n=l 

we can see that the four terms in Equation (3.8), which are composed of an even number of 
the helical PMF power spectra and proportional to or AbA'^, give the signals under the 
condition as 

3 

J2i^n + Xn) =eYen . (3.12) 

n=l 

These signals can arise from the parity-even CMB bispectrum as 

AX(^)(n)\ AX(^)(-n)\ , , 

and are produced even in the absence of the helical PMF, namely, Ajg = [28-31]. On the 
other hand, the other four terms in Equation (3.8), which are composed of an odd number 
of the helical PMF power spectra and proportional to A'^Aq or A^, contribute the signals 
under the condition as 

3 

5^(4+Xn) =odd . (3.14) 
n=l 

These signals are due to the parity-odd CMB bispectrum as 

{n^)^-{n^). 

and can not appear unless there exists the helical PMF. 

As discussed in Ref. [28], Equations (3.5) and (3.8) also have the one- loop structure in 
terms of multipoles, which is originated from the sexuplicate dependence on the Gaussian 
PMFs, and hence the numerical calculation of the CMB bispectrum requires a great deal of 
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time. To speed up the numerical computation, we can use the optimal formula generated 
from Equation (2.11). In the same manner as the formulation of the exact bispectrum (3.8) 
[28], we gain the approximate formula for the CMB bispectrum of the scalar and tensor 
modes as 



\n=l 



app 



R.y ln{Tiy/TB) 



L,L'L" 



ub + 3 * 3 



E (-1) 



-'L1L2L3 



L L' L" 
1 1 1 



rOl-l rOl-l 7-01-1 




v93-JVs r|Ai|0-|Ai| .|A2|0-|A2| ;-|A3|0-|A3| 



tiLiL ^i2L2L' 
3 



x8 / y^dy 
Jo 

y^jriZiZiZi) 



Hi-')' 



.n=l 



kldkn^(z^) 



27r2 



T'x:l{kn)jLAkny) 



(3.16) 



with 
-piSSS) 



X {Pi,(fc2)Pg,,„.,(5L',0 + 6l>,2) + P6(A;2)P[t,„./LM} V^'^ 
PBikl)V^LL,^^{SL,0 + Sl,2) + ^'6(fel)<U,./L,l} 



-^'3/3.2:3 



6l" 



X</2,X2'^^'.2 {^b(A.3)<^,,3,.3(5l",0 + 5l",2) - PB(^^3)^il.3,:.3'^^'M} 



+<A,./l,2 {^'B(/C2)I^it,.,(5L',0 + 5l',2) - P6(A^2)P[t,./LM} 
X 



{PB(/C3)2?g,,3,,.3(5L",0 + 5l",2) + P6(A^3)<^,,3,.3 5l'm} 
^(^^^) = 6^/3 [{PB{k,)vil^^^^{SL,o + 6l,2) - PB{kl)V^LL^^6L,l} 

X {PB{k,)V^^l^^^^{6L',0 + '^L',2) + ^'6(A;2)Pg,,,./LM} ^^^^ 



(3.17) 



-J^3/3,a:3 



+3 { (Pb(A:i)p£,^,,, (5l,o + <^L,2) + ^Mfci)<^A,./L,i) <^,.„./l',2 

(PB(A:2)P[t,.,(<^L',o + 5l',2) - ^Mfc2)I^g,,,../LM) } 



{PB(fc3)<,3,3-PM^3)<,3,3} 



2 1 



(3.18) 



^(^^^) = 108 [{Pb(^^i)<,,.,(5l,o + Sl,2) - ^'B(fci)<,,./L,i} 



V 



is /s, 2:3 



X (Pi.(A:2)Pr;,,,., 

PBikl)V^L!A,,^ih,0 + 5l,2) + ^'Mfcl)<^A,./L,l} 



XPI!).,.., (PMfc3)<,3,.3 - ^Mfe3)Pit 



^2/2,^2 ' 

35i,2 {PBik2)V^Ll 



,X3 



.,.2 - PB{k2)VZ,2,.2 



PB{k3)V^LL,., 



(o) I 



^L'2^L'\2 



(3.19) 
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^(^^^) = 648V3 



^2 I L3,£3,a;3 



Here, V^^\ ^ 



{o) 

Li,ei,xi 



(defined by Equation (3.10)) and 



(3.20) 



~2 + SLi,£i + SLi,ei+2)Sxi, 



(3.21) 



wliich lead to Li+ ii + xi = even and = odd, respectively, and the Kronecker delta functions 
arise from the selection rules of the Wigner symbols. Note that althought the vector modes 
have not been considered due to their negligible signals at large scales, we can formulate them 
in the same way. These functions and a selection rule as Yln=i = even ensure that the 
signals of the CMB bispectra satisfying Yl\=i{^n + a^n) = even are sourced from the terms 
proportional to and AbA^. This is consistent with the discussion in the exact formula 
(3.8). Likewise, we can see that the signals under X]n=i(^" •^") ~ '^^^ proportional to 
only A^Aq. This means that the contribution of A^, which appears in the exact formula 
(3.8), is negligible because an integral at around each pole in the bispectrum of the PMF 
anisotropic stresses is independent of the helical PMF as seen in Equation (2.12). 

Through numerical calculations, we confirmed that these optimal formulae reconstruct 
the shapes of the CMB bispectra based on exact formula (3.8) if a's are identical to the 
values for the non-helical case [28] as 



a 



0.3350 


{Zi 


= Z2 


= Z3 


= S) 




0.3473 


[Zi 


= Z2 


= s, 


Z3 = 


T) 


0.3212 


(Zi 


= s, 


Z2 = 


Z3 = 


T) 


0.2991 


(Zi 


= Z2 


= ^3 


= T) 





(3.22) 



3.2 Analysis 

Here, we show the numerical results of the CMB bispectra and signal-to-noise ratio. To calcu- 
late the CMB bispectra, we modified the Boltzmann Code for Anisotropies in the Microwave 
Background [42, 52] and used the Common Mathematical Library SLATEC [53]. 

Figure 1 presents the parity-even signals from Y^=i = even and parity-odd ones 
from Ylf^^i in = odd in the reduced intensity-intensity-intensity (///) bispectra of the three 
tensor (TTT), two tensor and one scalar {STT + TST + TTS), and one tensor and two scalar 
{SST + STS + TSS) modes defined by 



,{ZiZ2Zs) 

^III,ili2h 



l2ti) 



h h h 
mi 771-2 'IT'S 




(^n) 



(3.23) 
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Figure 1. Absolute values of the CMB /// bispectra generated from the TTT (red solid line), 
STT + TST + TTS (blue dashed one) and SST + STS + TSS (magenta dotted one) modes as the 
function in terms of £3 when the two multipoles arc fixed as (£i,i!2) = (100, 105). The curves in the 
left and right panels correspond to the parity-even and -odd bispectra arising from X]n=i = even 
and = odd, respectively. The PMF parameters are fixed to BiMpc = 4.7nG,SiMpc = 1.35nG,nB = 
?7B = —2.9 and r^/r^ ~ 10^^, and other parameters are identical to the mean values derived from the 
WMAP 7-yr data [54] 



where 




^^i^2^3 - ^^(^^ ^ 1) _ ^^(^^ ^ 1) _ ^^(^^ + 1^ V 4^ I -1 1 j • ^"^-^^^ 



This symbol is identical to I^^i^^.^^ when X]n=i ^™ ~ even [12, 55]. Here, we focus on the signals 
when multipoles have different values as (^1, ^2) = (100, 105) because the CMB bispectra from 
Yl\=i = odd are exactly zero due to their asymmetric nature in the case where ^1 = -^2 = ^3 
is satisfied [12]. Here, since we have fixed the strengths and spectral indices of the non-helical 
and helical PMF as Psik) = Psi^) is satisfied, it is natural prediction that the parity-even 
and -odd /// bispectra have almost same magnitudes in each mode. Actually, however, the 
parity-odd signals are smaller than the parity-even ones in each mode because the parity-odd 
signals are highly damped as three multipoles approach the similar values. The parity-even 
signals have same features as those in the non-helical case [28]: the TTT mode dominates 
for £3 < 100 due to the Integrated Sachs Wolfe effect [56] and the bispectra including the 
scalar modes gradually increase for £3 > 100 by the acoustic oscillation. In contrast, it 
is hard to observe these features in the parity-odd signals due to the damping effects at 
£^ £2 ~ ^3 ~ 100. 

Figure 2 shows the signal-to-noise ratios from the parity-even and -odd CMB /// bis- 
pectra of the TTT, STT + TST + TTS and SST + STS + TSS modes, and their total 
bispectrum, respectively ^, given by 

2 



S _ ^ , 



where A^^^j^ = 1 (fo^ ^1 7^ ^2 7^ ^3), 6 (for £1 = £2 = £3), 2 (otherwise) denotes a numerical 
factor caused by the cosmic variance and Cf'^ is the fiducial CMB power spectrum consistent 



We neglect the three scalar mode for its weak signals at < 500. 
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Figure 2. Noise- free signal-to-noise ratios from the parity-even (left panel) and -odd (right one) 
CMB III bispectra coming from J2n=i = even and = odd, respectively. The "total" line denotes 
S/N obtained from the total spectrum of the TTT, STT, TST, TTS, SST, STS and TSS modes, 
and the others correspond to S/N^s coming from each mode. Here, we fix the PMF parameters as 
-BiMpc = l-0nG,i3iMpc = 0.287nG, ns = = —2.9 and t^/tb = 10^''. The other parameters are 
identical to the mean values obtained from the WMAP-7yr data [54]. 



with the current observational data [54]. Here, we have neglected the noise other than the 
cosmic variance. We can find that both in the parity-even and odd bispectra, the TTT mode 
dominates over the signals. However, the scalar mode increasingly afltects at intermediate 
scales, where the tensor mode reaches a plateau, and will become important at smaller scales. 
Although we also have taken the PMF parameters keeping Pb(A;) = Pg(A;), the parity-odd 
signals are small compared with the parity-even ones. Nevertheless, these signals are precious 
information to probe the parity violation of PMFs because these should vanish if ^BiMpc = 
at variance with the parity-even signals. As discussed in the previous subsection, the CMB 
/// bispectra (xi = X2 = = 0) of the scalar and tensor modes from Yln=i = odd 
are asymptotically proportional to Bf^p^Bf-^p^[ln{T,^/TB)]^- Therefore, from the value of the 
total signal-to-noise ratio for £max = 500 described in Figure 2 as S/N = 2.236 x 10~^, we 
can understand that if ub = = -2.9, t^/tb = W^'^iW^) and ^iMpc^mpc > 2.7(4.5)nG, 
S/N exceeds unity and the existence of the helical PMF will be observationally supported. 



4 Summary and Discussion 

In this paper, we investigated the effects of the parity-violating helical PMF in the CMB 
bispectrum. At first, following the definition of the helical PMF in Ref. [23], we calcu- 
lated the bispectrum of the PMF anisotropic stresses involving both non-helical and helical 
PMFs. Then, through the application of the pole approximation [28], we found that the non- 
Gaussianity of the PMF anisotropic stresses closes to the local-type one if the non-helical 
and helical PMF power spectra have nearly scale-invariant shapes. In the same manner 
as our formalism [28, 29], we formulated the CMB bispectrum generated from not only 
non-helical but also helical PMFs. Our formulae showed that the finite signals arise from 
not only Yln=i ~ ^ven but also Yln=i ~ ^'^'^ intensity-intensity-intensity bis- 

pectrum. The latter signals, which is asymptotically proportional to B\y[^J3\y[^^, never 
appear unless the helical PMF exists and hence these are good observables to probe the 
parity-violating information of the PMF. Through the numerical computation of the CMB 
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intensity-intensity-intensity bispectra of the scalar and tensor modes and their signal-to-noise 
ratios, we found that the parity-odd signals from ^^=1 = odd are smaller than parity- 
even ones from Y^n=i ^" ~ even since the parity-odd signals are highly suppressed in the case 

where ~ £2 ~ ^3 due to the asymmetric nature of the CMB bispectrum. Nevertheless, the 

2/3 1/3 

computation of the signal-to-noise ratio provided a fact that if J^^f^^JS-^y^^^ > 2.7 — 4.5nG 
is satisfied, the parity-odd signals dominate the cosmic variance and it is possible to access 
the parity-violating nature of the PMF. Furthermore, the information of the polarizations is 
expected to improve this detect ability. 

This study with Ref. [12] gives a motivation to constrain the parity-violating non- 
Gaussianities from the observational data by using the signals under Yln=i ~ °dd in the 
CMB intensity-intensity-intensity bispectrum. In Ref. [12], we have shown that the parity- 
violating Weyl cubic term produces not the local-type non-Gaussianities but the equilateral- 
type ones in the graviton sector and they also induce the CMB intensity-intensity-intensity 
bispectrum for Yln=i = odd. Therefore, it will become important to constrain the magni- 
tude of each- type non-Gaussianity like f^^^, /nl"'' /nl'^"^ ™ parity-even case [54, 57], 
and differentiate between various non-Gaussian sources involving the parity violation. These 
considerations remain as future issues. 
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A Projection tensors 

In this section, we present some useful relations of the projection tensors based on Refs. [28, 
41]. These are utilized as tools for the formulation of the CMB bispectra (3.8) and (3.16). 

An arbitrary unit vector, a normalized divergenceless vector and transverse-traceless 
tensor can be defined by using the spin spherical harmonics as 

m 

6(±i)(k) = T5]C±in™(k) , (A.i) 

m 

4f^(k) = ^^4^^Hk)4^^^(k), 



with 




i {5m,i + 5^,-1) I , (A.2) 



V25, 



'm,0 
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(A.4) 



whose scalar products are given by 

I Att I An 

Then, the divergenceless vector and transverse-traceless tensor obey 

6W(k)ei^')(k) = 5A,-v (for A, A' = ±1) , 
e(f)(k) = fc„elf)(k) = 0, 
eif^k) = eif^(k) = eif)(-k), 
4? (k)eif (k) = 25,,_ V (for A, A' = ±2) . 

Using the divergenceless vector, a projection tensor in Equation (2.3) and a unit vector 
are expressed as 

p,,(k)= 4'^Hk)4"''^(k) , 

'^=±1 (A.5) 
^c = nMc6i+^^(k)4"'^(k). 

These relations lead to Equation (2.5) 

The projection tensors, which are expanded by the spin spherical harmonics as 



1 

— ( 
3 

■ 2 1 1 
M nia nif, 



= -2lll^' ^ y2M(k)arar 

Mmami, 

02"(k)sttf"'(k)+4ef"(k) 

Mma'mf, 

.(±2).,;^ _ ^(±2), 



oir(k)=eir(k) 

^ I M ma mf, 



2V3I°i\-i ^2y2M(k)a^ar ' ^ 1 1 



MrrLaTrn, 



decompose an arbitrary physical tensor such as the metric or the PMF anisotropic stress into 
the two scalar (Xisoj X '''''*)) two vector (x^^"^^) and two tensor (x^^^^) components: 



Xa6(k) = --Xiso(k)(5afe + x^"'(k)0:^^,^(k) 

+ Y x(')(k)0iJ)(k) + Y X^')(k)02)(k) . (A.7) 

A=±l A=±2 
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Here, = y^3/(87r) is given by Equation (3.7). Considering Equation (A. 4), we can 

derive the inverse formulae as 

x(°)(k) = ^Oi?(k)Xa.(k), 

x(±^Hk) = io(f)(k)xa.(k), (A.8) 

X(=^^)(k) = i0if)(k)xa.(k) . 

These are used in the calculation of the each- mode initial perturbation (3.3). 

Finally, we describe other useful expressions for computation of the CMB approximate 
bispectrum (3.16) as 

oi?(k)^^(k) = ^4'^)(k), 

4i^(k)4'^^(k) = ^^e^^(k)5^^_. , (A.9) 

4^Hk)4-^^(k) = -E^^rr E ■ 
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